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Abstract 

We report results of molecular dynamics simulations of a binary Lennard-Jones system at zero 
pressure in the undercooled liquid and glassy states. We first follow the evolution of diffusivity and 
dynamic heterogeneity with temperature and show their correlation. In a second step we follow 
the ageing of a quenched glass. As diffusivity decreases with ageing, heterogeneity increases. We 
conclude that the heterogeneity is a property of the inherent diffusion of the relaxed state. The 
variations with aging time can be explained by annealing of quenched defect structures. This 
annealing has the same decay constants for both diffusivity and heterogeneity of both components. 
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I. INTRODUCTION 



Diffusion in glasses and their melts is fundamentally different from diffusion in crystalline 
lattices. It has been studied in experiment intensively for many years. The increase in 
computer power in recent years now makes computer simulation studies near and slightly 
below the glass transition temperature possible. Combining the results from experiment and 
simulation one increasingly gains insight into the underlying atomic dynamics. Despite this 
effort, many questions remain and no general agreement on the nature of diffusion and its 
change at temperatures near the glass transition has been reached, even for simple densely 
packed glasses, such as binary metallic glasses. However, the research of recent years has 
greatly improved our understanding and several major issues have been settled; see the 

n 

recent review 

In a hot liquid, diffusion is by flow, whereas, in glass, well below the transition temper- 
ature, it will be mediated by hopping processes. Key question are the transition between 
the two regimes and also the nature of the hopping process. Are the jumps governed by a 
vacancy mechanism, similar to diffusion in the crystalline state, or by a mechanism inherent 
to the disordered structure? 

Glasses are thermodynamically not in equilibrium; their properties depend on the pro- 
duction history and one observes ageing. The diffusion coefficient of a glass which has been 
relaxed for a long time will be considerably lower than the one of an "as quenched" glass. 
Experimental evidence suggests that this could involve a change of diffusion mechanism. 

In the glassy state, diffusion can be described by an Arrhenius law. The activation 

n 

enthalpies are typically 1 to 3 eV jlj. Other than in crystals the pre-exponential factor 
varies widely from about 10~ 15 to 10 +13 m 2 s _1 l] 
often described by a Vogel-Fulcher-Tammann law 



I. In the undercooled liquid, diffusion is 

LL 



D{T) = D^ FT exp-H/k(T - T VFT ) (1.1) 

which describes a stronger than exponential decrease of diffusivity upon cooling toward 
T VFT < T g where T g is the glass transition temperature. As T — ► T VFT , a residual hopping 
diffusion eventually becomes essential and the diffusion crosses over to the Arrhenius law. 
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Mode coupling theory gives a different expression for diffusion in the undercooled liquid 

D(T) = Df CT (T -T c f . (1.2) 



The critical temperature T c , where viscous flow is arrested according to this expression, lies 
between the glass transition and melting temperatures, T VFT < T g < T c < T m . Again, 
hopping terms not included in the simplified expression will become dominant and prevent 
this freezing. The diffusion data alone can be fitted by both expressions. Neither expression 
allows a clear insight into the atomic process behind diffusion. 

From the pressure dependence of diffusion one finds, in general, apparent activation 
volumes much smaller than an atomic volume Q. For vacancy diffusion the activation 
volumes are of the order of the formation volume, i.e. of the order of the atomic volume. 
Low activation volumes hint at diffusion without formation of defects, i.e. diffusion by 
an inherent mechanism. A key to the nature of diffusion was found in measurements of 
the isotope effect of diffusion E. It is defined by the ratio of the diffusivities, D a , Dp of 
two isotopes with masses m a and mp, respectively E = {Da/Dp — 1) / {^J vn a jmp — j^. 
Because of the 1/^/m- dependence of the atomic velocities E is of order unity for diffusion 
via single vacancy jumps in densely packed lattices where essentially single atoms jump jf|. 
In contrast both in glasses Q and in undercooled metallic melts |8 : ] nearly vanishing isotope 
effects were observed. The vanishing isotope effect has a natural explanation if one assumes 
a collective diffusion mechanism both above and below the glass transition. The usage of 
the term collective follows the general usage in the field, see e.g. [1]. It indicates that groups 
of atoms jump together as opposed to the usual jumps into vacancies in lattice diffusion. 

The vanishing isotope effect is a property of relaxed glasses. Upon quench, defects can be 
frozen in. These enhance diffusion and can lead to an increased isotope effect which vanishes 
upon relaxation [jj. Quasi vacancies can be produced also under irradiation which again 
enhances diffusion, see |ic| . 

In molecular dynamics simulations one can follow the motion of atoms over periods of 
several ns and up to /is. Early simulations were restricted to simple model systems, such as 
Lennard- Jones [n| or soft sphere systems ^1 above the glass transition. Now it is possible 
to simulate real systems such as e.g. NiZr 13] or CuZr [l^. As far as qualitative properties 
of metallic glasses are concerned the results of the model systems and the real ones are in 
full agreement and can be used interchangeably. 

Extensive studies of different properties of an undercooled binary Lennard- Jones system 
(BLS) at constant volume showed good agreement with the predictions of MCT. For 

the pressure derivative this theory predicts a singularity at T c which should be modified by 
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hopping processes to a sharp cusp of the apparent activation volume at T c . This was observed 



in a simulation of a BLS at zero pressure 17[. Both above the cusp, in the undercooled 
liquid, and below, in the glass, activation volumes of around 0.4 atomic volumes were found. 
In the hot liquid the activation volume rises to values near the atomic volume, which is 
expected for binary collisions. 

The change of the diffusional isotope effect upon cooling was studied at zero pressure 

n Q 

for monatomic Lennard- Jones liquids |l8( and for BLS |l9j ]. In a hot liquid one has values 
E « 1. Upon cooling, E drops to below E = 0.1 already well above both T g and T c . This 
drop is mainly driven by densification but at different rates for the two components of the 
BLS. 

The small isotope effect can be understood from the elementary step of diffusion, the 
atomic hopping. In a soft sphere glass, at low temperatures, this atomic hopping was found 
to be highly collective. Chains of ten atoms and more move together in a jump. Whereas 
the total jump length, summed over all participating atoms, is of the order of a nearest 
neighbour distance, a single atom only moves a fraction thereof j^ol . 2l| . The same chain 
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like motion was observed in snap shots of the motion in the undercooled melt 
This collective motion is similar to the one of a group of people threading their way through 
a crowd. They will fill whatever space they find in front - small activation volume, and each 
member will follow instantaneously the one in front - collectivity. Obviously they will get 
on better as a file than in a broad front. Chain like motion was observed already in early 
simulations of melting in two dimensions j^ . 

Subsequent atomic jumps are strongly correlated with each other, not only in the case of 



return jumps 



2l|, and atoms which have just jumped have a strongly enhanced probability 



to jump again, leading to local bursts of activity interspersed with times of relative calm 
|2lL [26]. These correlations, which are typical of all glasses, not only metallic ones, are 
closely related to the so called "dynamic heterogeneity" of glasses and undercooled liquids. 
At any given time only a few atoms are active. The resulting mean square displacements 
strongly deviate from a Gaussian distribution for long times (non-Gaussianity) before it is 
finally restored by long range diffusion. 

In this paper we will present results for the temperature dependence of the dynamic 
heterogeneity and then study effects of ageing on both the diffusivity and the dynamic 
heterogeneity. 
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II. SIMULATION DETAILS 



The calculations are done for a BLJ system 

VijiR) = Aeij [((Tij/R) 12 - K-/.R) 6 + A i3 R + B i3 ] . (2.3) 

where the subscripts ij denote the two species. The potential cutoff was set at R c = 3a. 
For the parameters we took the commonly used values of Kob and Andersen |l5j: €aa — 
e = <jaa = u = I, €bb = 0.5, a BB = 0.88, €ab — 1-5 and oab = 0.8. To avoid spurious 
cutoff effects we introduce the parameters and Bij to ensure continuity _of_ the potential 



and its first derivative at the cutoff, similar to the shifted force potential [27[. All masses 
were set to rrij = 1. As usual, in the following, we will give all results in the reduced units 
of energy €aa, length oaa-, and atomic mass uia- To compare with real metallic glasses one 
can equate one time unit ((e/mcr 2 ) -1 / 2 ) roughly to f ps. 

The calculations were done with zero pressure and periodic boundary conditions. The 
time step was At = 0.005. Control runs with At = 0.0005 showed no significant deviation. 
A heat bath was simulated by comparing the temperature averaged over 20 time steps 
with the nominal temperature. At each time step 1% of the temperature difference was 
adjusted by random additions to the particle velocities. Apart from the very first steps of 
the ageing procedure the correction, after excursions of the temperature due to relaxations, 
did not exceed 10 -4 of the average velocity. This procedure assured that existing correlations 
between the motion of atoms were only minimally affected. The results did not change when 
the temperature adjustment was varied within reasonable limits. 

The ageing studies were done for samples quenched to T = 0.32 e/k, below the glass tran- 
sition temperature, defined from the kink in the volume versus temperature curve. Ageing 
leads densification. However for the parameters adopted the aged system was still suffi- 
ciently far from the density of the undercooled liquid to be considered glassy. We proceeded 
from the samples prepared in our previous work on the pressure dependence of the diffu- 
sivity where we had prepared 16 samples for each temperature, 8 with slightly positive and 
8 with slightly negative pressure Each sample consists of 5488 atoms in a ratio 4 : 1 
of A- and 5-atoms. The samples were then additionally aged for up to 32- 10 6 time steps. 
Measurements were done during runs with constant volume. In all runs pressure and energy 
were monitored to ensure stability of the configurations. The measured pressure was used 
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to interpolate to zero pressure. The diffusion constant was calculated from the asymptotic 
slope of the atomic mean square displacements. 

Fig. n shows the densities and diffusion constants for zero pressure. From the change in 
slope of the volume expansion we estimate the glass transition temperature as T g w 0.35e/fc. 
The diffusion constant can be fitted very well with the mode coupling expression, Eq. 11.21 
using a value T c = 0.36e//c for both species and 7 = 1.87 and 7 = 2.02, respectively. The 
two temperatures are very close to each other, T g w T c , but are much lower than the value 
T c = 0.435e//c, reported for simulations at constant density p = 1.2 This reflects the 
strong dependence of the glass transition on density or pressure. We find for zero pressure 
a density of p = 1.16 at T g . 

The diffusion coefficients depend in a rather intricate way both on temperature and 
atomic density. This makes a comparison of the present zero pressure values with the the 
previous isochoric ones [15j difficult. First the BLS becomes at p = unstable for T le/k 
whereas the high density used in stabilizes the BLS to T > 51e/k. As mentioned in the 
introduction the pressure derivative of D (activation volume) has a cusp at T c this implies 
that one would have to scale temperature as T/T c and subsequently scale with pressure 
which again is complicated by the strong temperature variation of the activation volume 
which reflects the transition from a thin liquid dominated by binary collisions through the 
viscous undercooled liquid to the solid like glass jl?] ]. Furthermore it has been shown that 
both components are affected differently by density. Density is a strong driving force towards 
cooperative motion. However, in a binary liquid there is no longer a simple proportionality 
jlfl| . We have, therefore, not attempted to scale our values over the whole temperature 
range to the ones of Q Doing a rough scaling just above T c we find qualitative agreement. 

Notabene, the present values for the diffusivity below T g , in the glass, are somewhat lower 
than the ones published previously (17J which is an effect of the additional ageing as we will 
discuss further down. 



III. DYNAMIC HETEROGENEITY 



In isotropic diffusion the atomic displacements are Gaussian distributed. In undercooled 
liquids and in glasses Gaussianity is violated over long time scales. This non-Gaussianity 
indicates different mobilities of different atoms over long time scales, the so called dynamic 
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heterogeneity. This effect is quantified by the non-Gaussianity parameter (28 ] 



3 < Ar 4 (t) > 1 . . 

where < ... > denotes time averaging, Ar 2 (t) is the mean square displacement and Ar (t) 
the mean quartic displacement. This parameter is defined so that it is equal to zero when the 
atomic motion is homogeneous. Experimentally it can be obtained from the g-dependence 



of the Debye- Waller factor 
simulations of liquids, e. g. 12, 



It has been calcu 
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ated in numerous molecular dynamics 
331 ] . There are three time regimes of 



0.2(f)- First it increases from a 2 (t = 0) = on a ps time scale due to heterogeneities of 
the atomic vibrations. The maximal value of 02 (t) in this vibrational regime is small, less 
than 0.2. Depending on temperature this is followed by a strong increase during the time of 
the so called /^-relaxation. At times, corresponding to the a-relaxation time, 02 (t) reaches 
a maximum and drops for t — > 00 to 012 f = 00) = 0. The strong increase seems to follow a 
Vt-law, independent of the material. This Vt-law has been shown to be a direct consequence 
of the collectivity of the diffusional jumps and the correlation between subsequent jumps 
j^ . This general behavior holds both below and above T g . 

In Fig. El we show the maximal value of the non-Gaussianity as function of inverse tem- 
perature. In the hot liquid a2(t) mauX is 0.12 and 0.14 for the two components, respectively. 
This corresponds more or less to the vibrational contribution with very little addition from 
jump processes. The hot liquid is, as expected, nearly homogeneous with respect to diffu- 
sion. Cooling down, we observe a marked increase of ce2(t) mauX i n the undercooled melt which 
accelerates approaching T c . At T c the maximal non-Gaussianity is already 2.5 and 4 for the 
two components, respectively. In the glass, just about 10% below T c , these values have dou- 
bled and reach a>2(t) max ~ 10 for the smaller component. This value is still an underestimate 
due to ageing effects, see next section. The non-Gaussianity is strongly pressure dependent. 
We define a pressure coefficient as 

Pa(T) = -L- - a2 Z ,Pl \~ a2 Z ,P2 l (3-5) 
P1-P2 a 2 (T,pi) + a 2 (T,p 2 ) 

with pi and P2 two different applied external pressures. We find for T = 0.32 f3 a (T) = 
0.27 and 0.35e/a 3 for the two components, respectively. There is no drastic effect near T c 
(j3 a (T) = 0.27 and 0.35e/(x 3 ). 

The increase of the maximal value is concomitant with an increase of the time this value 
is reached, t max (T) Fig. 01 In the hot liquid t max (T) is of the order of vibrational times and 



increases by four orders of magnitude upon cooling to 0.9T C . This slowing down reflects the 
general slowing of the diffusional dynamics as shown in Fig. ^ To check the correspondence 
between the diffusion time and t max we calculate the dimensionless quantity 

C D _ NG (T) = D(T)p(T) 2 /H max (3.6) 

where p(T) is the density of the system at zero pressure. In defining Cd-ng(T) it is assumed 
that the non-Gaussianity is mainly caused by the same atomic motion as diffusion which is 
indicated by the rise of t max (T) above the vibrational times upon undercooling. Cd-ngCO 
should then become independent of temperature. Indeed the large variation of t max (T) by 
orders of magnitude nearly vanishes. We find, for both components, C*d-ng(^ = 0.56) ~ 
0.015, at the onset of undercooling, dropping to Cd-ng(^ = 0.32) ps 0.004 just below T c . 
This correlation between t max (T) and diffusivity is related to the one with the onset of the 
a-relaxation, reported earlier [lfij j. The drop of Cd-ng(^) is probably partially a result 
of the increasing separation of relaxation and diffusion, the first one being less sensitive 
to eventual return events than the long range diffusion.. Another contribution could be a 
change of the shape of the non-Gaussianity versus time curves. 

As mentioned in the introduction, there is a general -\/(t)-law governing the increase of 
ct2(t) above its vibrational value 2^. From this, one could assume that a^t^Vax/ \/t mSiX (T) 
should be approximately constant in the undercooled and glassy regimes. This holds fairly 
well in the undercooled melt above T c but breaks down below. Part of this might be due to 
ageing but inspection of the actual «2(t)-curves shows that it is mainly due to an increased 
flattening near the maximum which separates t max (T) from the y^(t)-increase. 



IV. AGEING 



We have seen that there is a close connection between diffusivity and dynamic hetero- 
geneity. It is well known that glasses are not in thermodynamic equilibrium and are subject 
to ageing. In experiment one sees upon rapid quenching a drop of the diffusivity on a time 
scale of hours |24|. Ageing, it is generally assumed, leads to a more "ideal" glass. The 
"excess volume" drops which indicates annealing of defects - whatever they are. This poses 
the question of the relationship of dynamic heterogeneity and idealness. If the dynamic het- 
erogeneity is an inherent property of the glassy state one expects it to increase with ageing 
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as the diffusivity decreases. 

A simulation over real times of hours is of course impossible. To circumvent this, ageing 
effects have been studied by instantaneous quenches from high temperatures, thus producing 



large effects |35j,|36(. We adopted a softer procedure: samples equilibrated at T = 0.34e//c 
were quenched at a rate of Q = 10~ 4 to T = 0.32e/k. The quench amounted to about 5% 
of T c . The quenched samples were then aged and the diffusion constant was determined 
in subsequent intervals. To determine the diffusion constant at time t we calculated the 
average mean square displacement in the interval [t — 5/2, t + 5/2]. The diffusion coefficient 
was then obtained by the standard procedure from the long time slope in that interval. 

Fig. H] shows a drop of both diffusivities by about an order of magnitude over the ageing 
interval of 80000 time units. The rapid initial drop of log D(t) is followed by a seemingly 
linear part, dashed lines. Such a behavior is consistent with the assumtion that the diffusivity 
is the sum of two terms 

D(t) = Anh + AlcfCdcf (*) (4.7) 

where Anh is a time independent inherent diffusion coefficient, Aef is a defect mediated 
diffusion constant and c^it) is the defect concentration. If the defects are slowly annealed, 
with some rate «def during the ageing at constant temperature, we get 

D(t) = Anh + AefCdcf (0)e- ad < (4.8) 

Assuming one type of defect, the same decay constant a should apply to both components 
whereas the combination AefCdef(O) can vary between them. Eq. 14.81 gives, appart from the 
short time behaviour, an excellent fit of the data of Fig. HJ dashed line. The short time 
behaviour cannot be expected to be reproduced by such a simple model of independent 
defects. We find a decay constant ad e f(T = 0.32) = 4.15 • 10~ 5 , and for the combination 
AefCdef(O) we find values of 1.45 • 10 -6 and 5.84 • 10 -6 for the large and small component 
respectively. Eq. 14.81 allows us to extrapolate to the inherent diffusivity which should be 
reached after long time aging. These values are shown by the open symbols in Fig. ^ 

Soon after the quench, the first values of the diffusion coefficients equal within 10% 
the ones before the quench, at the higher temperature. There are two effects which seem 
to cancel each other. On the one hand, due to densification and temperature reduction, 
diffusion should slow down. On the other hand, the quench drives the system further away 
from equilibrium which 



accelerates diffusion processes. It should be noted that diffusion in this short time interval 
is not necessarily long range. 

During the ageing the glass is densified. The volume reduction per atom is AQ pa 2- 10 -3 n 
where Q stands for the average atomic volume at that temperature. There is no direct 
proportionality with the drop in diffusivity. By the time the rapid drop finishes (t pa 30000) 
only about 20% of the excess volume, AQ, is gone. The raised diffusivity can, therefore, not 
be attributed to simple quasi- vacancies but to more intricate "defects" . 

The slowing down of diffusion is accompanied by a monotonic increase of the non- 
Gaussianity, Eq. 13 A\ by a factor of two over the ageing period, Fig. EJ Immediately after the 
quench the values of c^Wmax are considerably below the ones of the starting temperature 
T = 0.34e/A;. This is again consistent with the above picture that the quench produces some 
"defects" which are annealed during aging. This is described in analogy to Eq. 14.81 bv 

a 2 (t) = 4 nh + a fc def (0)e-^'. (4.9) 

Assuming that the slowing down of diffusion and the increase of heterogeneity are caused 
by the same mechanism, we take for the decay constant the value obtained from the diffusiv- 
ity. The resulting fit is shown by the dashed lines in Fig. [5J Considering the obvious scatter 
the fit is again excellent. The resulting values are for the inherent value a^ h = 6.2 and 10.0 
and for the defect part a2 Cf Cdef(0) = —4.0 and —5.1 for the two components, respectively. 

The close connection between collective jumps and heterogeneity gives a hint of the 
possible nature of what we loosely call defects. Defect here does not mean a simple structure 
such a vacancy but more likely a local center of strain. These local strains can be relaxed 
by a less collective motion than the one in the relaxed glass. This was found earlier by 
racer experiments experiments which investigated the effect of ageing on the isotope effect 
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This additional motion apparently is more homogeneous than the one inherent to the 
relaxed glass. At the early stages or after a large rapid quench, collectivity might be fully 
destroyed for some jumps. In a simulation of a BLJ glass at constant density single particle 
jumps have been observed after a rapid quench [37 1. 

Finally, checking the correlation, Eq. 13.61 we find a drop by a factor of 2-3 during ageing, 
similar to the one found in the temperature dependence. We interpret this as an indication 
that the excess diffusivity produced the quench shows less decoupling between diffusion and 
relaxation than the inherent one. 
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V. SUMMARY 



We investigate by molecular dynamics the relation between diffusivity and dynamic het- 
erogeneity We use a binary Lennard- Jones like system as simple model of a metallic glass, 
respectively melt. This model system has been extensively studied, mostly at constant vol- 
ume. In this study the volume is varied to have zero pressure conditions. The diffusion 
coefficients in the melt are in accordance with mode coupling theory with a single critical 
temperature T c . In the undercooled melt and, even more pronounced in the glass, diffusion 
is strongly heterogeneous over long times. The distribution of atomic displacements deviates 
from the Gaussian distribution found for random walks. 

The non-Gaussianity parameter increases, from its small starting value due to vibrations, 
initially with a y/i-law before long range diffusion finally restores homogeneity. The time of 
the maximum non-Gaussianity is given approximately by the diffusion time. However, the 
correlation factor between the two times decreases systematically with reducing temperature. 
This may indicate a decoupling between relaxations and long range diffusion as has been 
observed in studies of viscosity versus diffusivity. Below T g the A/t-increase does no longer 
determine the maximal value of the non-Gaussianity since the maximum rapidly flattens. 

After a quench below T c we find a drop of the diffusion coefficients by an order of mag- 
nitude with ageing at constant temperature. This can be explained in terms of defects 
which are produced in the quench and are subsequently annealed. These defects lead to 
faster diffusion and lower heterogeneity. In this simple system, defect probably means some 
center of local strain which allows for a less collective motion. In a more complicated sys- 
tem "defects" could, e.g., be different molecular conformations This general picture is 
supported by the fact that in the present binary system the aging of both, diffusivities and 
heterogeneity, of both components can be described by a single decay constant. This can 
be used to extrapolate to the inherent diffusion coefficients of the ideal glass at the given 
temperature At the temperatures accessible to simulation this ideal glass would actually still 
be an undercooled liquid, alas with T < T c . The heterogeneity increases with ageing. We 
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conclude that it is an inherent property of the relaxed glass which is suppressed by defects. 
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Fig. 1 Diffusion constants (majority A-atoms, diamonds, and minority 5-atoms, cir- 
cles) and density (triangles) at zero pressure against inverse temperature (all in 
reduced units). The open symbols indicate the extrapolated inherent diffusion 
coefficients according to Eq. 14.81 The dashed lines show the fit with the mode 
coupling expression for atoms A and irrespectively. 

Fig. 2 Maximum value of the non-Gaussianity (majority A-atoms, diamonds, and mi- 
nority S-atoms, circles) at zero pressure against inverse temperature. The dotted 
line indicates T c . The full lines are guides to the eye only. 

Fig. 3 Time, t max , elapsed till the maximum of the non-Gaussianity is reached (major- 
ity A-atoms, diamonds, and minority 5-atoms, circles) at zero pressure against 
inverse temperature. The dotted line indicates T c . 

Fig. 4 Diffusion coefficient as function of ageing time (majority A-atoms, diamonds, and 
minority i?-atoms, circles). The dashed lines show the fit with the exponential 
annealing expression, Eq. 



Fig. 5 Maximal value of non-Gaussianity as function of ageing time (majority A-atoms, 
diamonds, and minority 5-atoms, circles). The dashed lines show the fit with 
the exponential annealing expression, Eq. 
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FIG. 1: Diffusion constants (majority A-atoms, diamonds, and minority 5-atoms, circles) and den- 
sity (triangles) at zero pressure against inverse temperature (all in reduced units). The open sym- 
bols indicate the extrapolated inherent diffusion coefficients according to Eq. 14.81 The dashed and 
dash-dotted lines show the fit with the mode coupling expression for atoms A and B, respectively. 
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FIG. 2: Maximum value of the non-Gaussianity (majority ^4-atoms, diamonds, and minority B- 
atoms, circles) at zero pressure against inverse temperature. The dotted line indicates T c . The full 
lines are guides to the eye only. 




FIG. 3: Time, t m ax, elapsed till the maximum of the non-Gaussianity is reached (majority A-atoms, 
diamonds, and minority S-atoms, circles) at zero pressure against inverse temperature. The dotted 
line indicates T c . 
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FIG. 4: Diffusion coefficient as function of ageing time (majority ^4-atoms, diamonds, and minority 
.B-atoms, circles), he lines are guides to the eye only. The dashed lines show the fit with the 
exponential annealing expression, Eq. I4.8I 
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FIG. 5: Maximal value of non-Gaussianity as function of ageing time (majority ^4-atoms, diamonds, 
and minority S-atoms, circles). The dashed lines show the fit with the exponential annealing 
expression, Eq. I4.9I 
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